from constants import *
import numpy as np
xx = 300
atom_wt = 197
rho = 19.3
density = rho * na / atom_wt * 1e6
print "density = ", density
print "ngold = ", ngold
lens = 8e-6
dx = lens / xx
dt = dx / c_light * 0.3
denspc = 1e13
dens = denspc / dx / dx
dens = 1e20
print "density = ", dens
omega = np.sqrt(dens * eV**2 / me / epsilon_0)
freq = omega / 2.0 / np.pi
period = 1.0 / freq
print "period = ", period
cycle = period / dt
print "need cycle = ", cycle
